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Abstract 

For fields that vary slowly on the scale of the lightest mass the logarithm of the 
vacuum functional can be expanded as a sum of local functionals, however this 
does not satisfy the obvious form of the Schrodinger equation. For if 4 theory we 
construct the appropriate equation that this expansion does satisfy. This reduces the 
eigenvalue problem for the Hamiltonian to a set of algebraic equations. We suggest 
two approaches to their solution. The first is equivalent to the usual semi-classical 
expansion whilst the other is a new scheme that may also be applied to theories 
that are classically massless but in which mass is generated quantum mechanically. 



For fields that vary slowly on the scale of the lightest mass the logarithm of the 
vacuum functional of a quantum field theory reduces to a sum of local quantities. The 
purpose of this letter is to construct the Schrodinger equation that acts directly on this 
local expansion. This reduces the eigenvalue problem for the Hamiltonian to a set of 
algebraic equations for which we suggest two methods of solution. One of them is the 
familiar semi-classical expansion, the other is a new approach that would apply even when 
mass is generated quantum mechanically in a theory that is classically massless. 

To begin with, consider a free massive scalar field theory in D+l dimensions. In the 
Schrodinger representation the Hamiltonian is H = — ~A + / d D x | (Vy? • Vy? + m 2 (p 2 ), 
(see | l | for a review of the Schrodinger representation in field theory). The unregu- 
lated Laplacian, A, is J dxd~ 2 /d~ip 2 (x), and the vacuum functional ty[<p] = (f\0) = exp — 
/ d D ^ipTip/2 = expWo with T = +V— V 2 + m 2 . The corresponding energy eigenvalue, 
E, is proportional to the functional trace of T and so diverges. The Laplacian needs to 
be regulated so introduce a momentum cut-off, 1/s, by defining 

d D p 5 2 



i/s (2tt) d 5ip(p)5ip(—p) 



where y?(p) = / <i D xy?(x) exp — ip ■ x. The vacuum energy density S = E/V is now 
well-defined and diverges as the cut-off is removed 



£ = h A < W ° = \ U„ (tpV / ^^ ~ {D + l (D+1)l , (2) 

Where k is the area of the unit sphere in D dimensions divided by 2{2n) D . 

If ip varies slowly in space on the scale of 1/m then Wq [p] simplifies to a local expression 
obtained by expanding T as m — V 2 /(2m) — (V 2 ) 2 /(8m 3 ).. to give 

Similar local expansions may be performed for the wave-functionals given by describ- 
ing fields coupled to gauge-potentials. A local vacuum-functional for Yang-Mills was first 
discussed in || where it leads heuristically to an area law for the Wilson loop via a kind 
of dimensional reduction. Quark confinement from dimensional reduction was discussed 
in p. Dimensional reduction resulting from the local nature of the Schrodinger vacuum 
functional has been generalised by Horiguchi et al to the case of the Wheeler-De Witt 
equation to look for a new phase in quantum gravity beyond the Planck scale M. We 
expect that the vacuum functional of the interacting scalar theory can similarly be ex- 
panded as a local quantity for slowly varying if. The vacuum functional satisfies the 
Schrodinger equation, the existence of which was shown by Symanzik,[^]. In principle this 
determines the coeflicents of the expansion. However the local expansion does not satisfy 
the obvious form of this equation because expanding in local quantities does not commute 
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with removing the cut-off, s — > 0, even for the free theory. To see this apply the regulated 
Laplacian A s directly to the local expansion (|3|) to obtain the energy density as 

J-AW l ° c =f -f^(™-^-^L)=Y ^ (4) 

2V 3 Jp><i/ S (2ir) D \2 Am 16m 3 "; ^ {m 2 s) n+D / 2 K> 

(where a n = fcm D+1 T(|)/(r(3/2 — n)Y{n+ l)(D + 2n))). This expression appears to have 
divergences of increasing order as s —>■ unlike @ which correctly gives the behaviour 
of the vacuum energy as the cut-off is removed. This is because our local expansion Wq 
is only valid for slowly varying ip, that is for </?(p) with support in p 2 < m 2 , hence @ is 
only valid for sm 2 > 1, i.e. for large s. To obtain information about the vacuum energy 
for small s from this expansion we must re-sum the series. First define the continuation 
of the vacuum energy to complex s by 

By choosing the square-roots to have cuts on the negative real axis this is analytic through- 
out the complex s-plane with the negative real axis removed. For |s|m 2 > 1 it has a large 
s expansion which is just (£§). Let C be a key-hole shaped contour running just under the 
negative real axis up to s = —s with s > 1/m 2 , around the circle of radius s centred 
on the origin and then back to s = — oo running just above the negative real axis. The 
integral 

J(A) = -i 7 / d Ms)e^ (6) 



2m Jc s 

may be evaluated using the large s expansion (|) as 



oo ^ / ^ v n+D/2 



n=0 



T(n + 1 + D/2) \m 2 / 



We can also evaluate the integral by collapsing the contour C until it just surrounds 
the negative real axis. By taking A to be real, positive and very large the integrand will 
be exponentially suppressed away from the vicinity of the origin so that the integral is 
determined by S(s) for small s. In fact /(A) ~ \( D + 1 )/ 2 Y^=q &n^~ n as A — > oo with 
constant &i up to exponentially suppressed terms, so that /(A) provides a regularisation 
of the vacuum energy (H). If instead we had been dealing with a function S(s) that 
was finite in the vicinity of the origin with poles and cuts only on the negative real axis 
we would have obtained S(0) up to exponentially suppressed terms. In an interacting 
theory we would only be able to compute a finite number of terms of these expansions. 
Now (0) is an alternating series so truncating it at n = N gives an error less than 

I r(Af+2+ 1 D/2)-m 2JV + 2 +- p I wn i cn f° r large iV behaves as (e\/ m 2 Ny N+l+D l 2 ^ / N . Since we want 
to take A large set A = iV/i 2 , then the truncation error goes to zero with large N provided 
the iV-independent mass-scale \i is smaller than the particle mass m. This shows that we 
can extract information about the high momentum cut-off theory by working to a finite 
order, N, with the local expansion of the vacuum functional valid for slowly varying fields. 
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We now turn to the interacting theory. Symanzik || showed that for D = 3 the 
Schrodinger functional is finite when the cut-off is removed provided that in addition 
to the usual renormalisation procedure there is a further field renormalisation to take 
account of divergences associated with the boundary. This was taken further in || . The 
appropriate renormalisation constants are computable in perturbation theory, but since 
the theory is not asymptotically free these are not reliable so instead we will work with the 
1+1 dimensional theory which is super- renormalisable. In the usual treatment without 
boundaries the only divergent Feynman diagrams for (p 4 theory in 1+1 dimensions are 
those in which both ends of a propagator are contracted at the same point. These may 
be eliminated by normal ordering the Hamiltonian with respect to a perturbative vacuum 
to obtain a regular operator that does not depend on a cut-off: 



H 



dx : 



\^{xf + ^\xf + M^{xf) + ^ 4 



Here M r and g are finite, and the normal-ordering is with respect to oscillators with mass 
M r . There are no further divergences in the Schrodinger functional due to the boundaries 
so in this case there is no further field renormalisation. (This is thought to apply also in 



Yang-Mills theories [lOfl). Now define an operator H s constructed from momentum cut 
off fields 



dx 



i(fr 2 + ^ + M 2 ( S )^) + ^ 



where 



and the kernel 



<p s (x) = / dyQ s (x,y)<p(y), n s (x) = / dyg s (x,y)n(y) 



(9) 



(10) 



IF 



lp 2 <i/s 2tt 

implements a momentum cut-off. Writing this in terms of normal-ordered operators en- 
ables us to express H as lim s | H s provided that 



and 



s = -f dp .^Tvr \ M2 ^- M r ) | 9 ([dp I 

Cs 2 J p2 <i/s 2vr I VP + m r + 2 / 2 + M 2 J + 32 [J 2vr y/p* + M\ 



dp 1 



Ml = M 2 (s) + 



4 JjKX/a 27T ^J p 2 + M 2 

so that the Schrodinger equation for the vacuum is just lim s j_ o ((^|/f s |0) = E((p\0) or 



lim 

40 



^A s + J dx Q ( V ' s (xf + M\s)^ s {xf) + |^ s (x) 4 - £{s))) * 



lim ( — A s + V S )V = lim H s ^ = 

s|0 2 slO 



(12) 



(13) 



(14) 



If we evaluate this expression for a ip of that has no Fourier modes with momentum greater 
than 1/s then we can replace ip s in this expression by ip itself. 
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We will now show that the limit of small s may be taken by re-summing the large 
s behaviour, just as for the vacuum energy of the free theory, enabling us to compute 
the coefficents an expansion of the logarithm of the vacuum functional in terms of local 
quantities, and the eigenvalues of the Hamiltonian. We begin by showing that A s \l/k/>] 
extends to an analytic function in the complex s-plane with the negative real axis removed. 
Consider (tp\7r(x, 0) tt(x', 0)|0), written in terms of a Euclidean functional integral. It will 
be helpful to make the ip dependence of (<p\ explicit. To do this we can use a bra (D\ 
which is annihilated by (p. (D stands for Dirichlet). Then 

{<p\ = (D\e i f dxv( - x) *( x) , (15) 

so that 

((p\ = i((p\Tt(x), (16) 



8tp(x 

and using the canonical commutation relations ((p\<p = <p(<p\- Now in Minkowski space 

T (n(x, t) n(x', 0) = (0(x, t) 0(x>, 0) - iS{x - x')5(t - t') (17) 

so that if we use the standard relationship between functional integrals and time-ordered 
expectation values we can express ((p\e~ lt ° H \D) as the functional integral 

V(p e iS l l p\+f dx ( i( P l P+ A <p 2 ) (ig) 

where \D) implies Dirichlet boundary conditions (p = at times and to- A is a regular- 
isation of the 6(0) that arises from (0). Similarly the time ordered Green's function 

T(ip\e- itofl 0(x,t)0(x',t')\D) = [ V(pe iS[ ®+$ dx ^ +K * 2) (p(xit)<p(x',t'), (19) 



Again making use of fllTD, rotating to Euclidean space and taking the limit t — > oo gives 
for > t > f 



(cp\7r(x,t)7t(x',t')\0) = 

- J Vfa-SBm+fMvv+Av 1 ) fa X} t ){p( x >, t ') - 6(x - x')6(t - £')) (20) 

where the integration variable (p is now defined on the Euclidean half-plane, t < 0, with 
the boundary condition <p(x,0) = 0. To obtain ((p\7t(x, 0) n(x', 0) 1 0) we can either take 
the limit as t, t' f keeping t > t', in which case the 6- functions will not contribute, or we 
can take t = t' = having first cancelled the 5-functions against a divergence coming from 
ip(x , t)ip(x f , t') . A similar argument applies to the A<p 2 terms. Whichever approach we 
take affects only the ^-independent and / dx <p 2 contributions to the connected diagrams 
in A s \l/. For example when we integrate against Q s to get A s \l/ the delta-functions give 
a term proportional to \P / 'y/s which contributes only to the l/\/s part of the vacuum 
energy in (114]). Apart from this the 5-functions have no effect, and henceforth we ignore 
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them. Setting t — t' — in the remaining term, rotate co-ordinates so that the two points 
(x, 0) and (x', 0) differ by a Euclidean time r = \x — x'\ to obtain 

J V£ r e- SErm+ I dt ^&(0,T)&.(0,0) (21) 

where now the integration variable <p r is defined on the half-plane x > 0. This can be 
interpreted as the time-ordered vacuum expectation value of fields that evolve in Euclidean 
time with a Hamiltonian H defined on the half-line x > 

T E (0r\ef dt * t) '? m <p'(0, T)<p'(0, 0)\0 r ) (22) 

The vacuum |0 r ) corresponds to the rotated Hamiltonian. Expanding the exponential 
and making the Euclidean time dependence explicit reduces this to a sum of terms of the 
form 

/ dt n dt n -i.. dt p dtp-!.. dt q dt q -!.. dtiip(t n )..ip(ti) 
Jt Jt Jt Jo Jo J—oo J—oo 

(O.^'e-^^-^^V^e-^^^V'e-^^^^V^-e-^^'e^-V^e-^^^VlO,) (23) 
The integrals may be computed if we first Fourier decompose <p(ti) as 

(p(U) = (27T)- 1 j dk t p(k)exp (ihU). (24) 

This gives a delta-function conserving the total momentum h and for t{ with i > p this 
gives insertions of (H — iJ2 h)" 1 whereas for % < q it gives insertions of (H + iJ2 h)' 1 . 
The remaining t-integrals may be done by inserting a basis of eigenstates of H between 
each operator which leads to a sum of products of energy denominators of the form 
(Ei — E2 — iJ2 multiplied by exponentials of r of the form exp — (E\ — iJ2 kj), so 
that the r dependence may be made explicit as a sum of integrals over the spectrum of 
H 

j dE dk 1 ..dk n (p(k 1 )..(p(k n )SQ2 hi) (po + Pie lkT + p 2 e* (fcl+fc2)T + ..) e' Er (25) 
Integrating this against the kernel Q s gives A s ty as a sum of terms like 

/ * / dEdk 1 ..dk n q>(k 1 ).rtk n )6(£,k i ) E ^^"vn ^ 

jp 2 <s-! 71 J ^ ^ E — i[p + 2^ ki) 

We will make the dependence on s more explicit by setting p = qj \fs. Also we will scale 
the function if with s by evaluating this for ip of the form (fi(x) = f(x/y/s) = <fi s (x) so 
that <p(k) = f(ky/s)y/s, and we take f(k) to vanish outside \k\ < k « 1. This enables 
us to scale the ki integrals to obtain 

J qa<i I / dE ikMkM^ k.) E Jf^- + - |L (27) 

The functions pj acquire a dependence on s via the energy denominators that can be 
written after scaleing of ki as l/((Ei — E 2 ) — i(J2k)/\/s). Up to this point s has been 
real and positive but the expressions we obtain for (A s , I')[0 s ] can now be continued as 
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analytic functions to the whole complex s-plane excluding the negative real axis, (at least 
order by order in the expansion in powers of (f) s that we used to treat the exponential in 
(p2p). Similarly ((H s — E)^>)[4> s ] extends to an analytic function of s with singularities 
only on the negative real axis, so its finite value at s = can be obtained from an integral 
over a contour on which \s\ is arbitrarily large. Thus the Schrodinger equation becomes 

lim -L / -e Xs ((H s - E)*)[<f> s ) = (28) 



We take C as before. Setting \1> = e gives 

fei/cT eA ' B i A - w + 1 dxdy5 ' ix ' y) ^m) - v -- e ) w = (29) 

where we have divided through by *&[(p s ] which has a finite limit as s J, as can be checked 
in perturbation theory. Since k « 1 we can replace Q s (x,y) in (|29|) by a delta-function. 
For a general W[tp] we would expect to be able to expand in powers of if to give 

°° r dk dk 

W[(p] = £ J ^..^m)--mn)KY,h)Tn{kl,.,k 2 n). (30) 

An expansion in terms of local quantities would result from expanding T n in positive 
powers of momenta for small momenta. Applying A s to W[ip] gives 

A s W[ip] = 

/, i X ) [ l^--^^^(h)-^(hn)5(J2 k i) T n(P^ -P,h,-,k 2 n)- (31) 

Jp 2 <\ 27T J In Z71 

For large s the momentum p is small so we can expand T n (p, —p, fc 3 , .., /c 2n ) in positive 
powers of p. If we evaluate A S W for if = (fi s then we can also expand in positive powers 
of the other momenta hi, again obtaining a local expression, since the Fourier transform 
of 4> s has only low momentum modes for large enough s. This coincides with what we 
would have obtained by applying A s directly to the local expansion of W followed by 
substituting (fi s for if. This leads to a set of algebraic equations which re-formulate the 
eigenvalue problem for the Hamiltonian. Concretely, we study 

hi Ic " S S " eA * ("5 + I dX (¥) 2 I ~ V '~ Ex \ W = (32) 



for large A taking W\ to be local. 

Integration by parts gives linear relations between local expressions so we pick as a 
linearly independent basis expressions of the form / dxip v ° (ip') Vl ..((p^) Vn where v n , the 
power of the highest derivative in the expression, must be at least two. So we take as our 
local expansion of W\ 



J dx(^p 2 (b + c ip 2 ) + <^ /2 (&i + ci^ 2 + d^' 2 ) 

oo 

+ £ ip {n)2 (b n + c n ip 2 + d n ip' 2 + e n ipip" + ..) 

2 

+V>\foV 2 + fiV' 2 + fiV" 2 + ••) + ■• ) (33) 
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The coefficents &j, q, .. depend on A through (|32|) but since \1/ is finite they are themselves 



finite as A — ► oo. Substituting this into our Schrodinger equation (R3) yields 



^, , 8M , 166 2 A 2 , , 26 nV ^FA n 

r(n + 3/2)(2n + 1) 



( 2£(A) + 2£ + ( Ab + ^ + + • • + WZ o I ^ + - 




,2 



2ciA c 2 A" 



2 / 2 I6M1 / 4c 2 A d 2 A 2 e 2 A 2 
+y? z — H =- H 2ci + 2diA H h 



7T \/7T 7T \ 3 3 



2 /2 1 

7= + 



/326 ci + 966ic 12/jVA 



\ V7T 7T / 

„ 2 /32& & 2 + 166? VA / e 2 A \\ 

+ ^ 1 3^ + v( C2 + ir-- + "JJ 

+.. 

/ (W)\2 / ^brrfin—m 

a/A /2c n 2(d fl — e„)A \\ 

j ^r(n + l/2) 7T \n\ 3(n + l)! 7/ 

+... J = (34) 

where ip = <p s for s = 1/A Requiring the coefficent of each linearly independent function of 
ip to vanish gives an infinite set of algebraic equations. These may be solved in the usual 
semi-classical approach by first ignoring the power series in A so that to leading order 
6 = vr 1 / 4 M/2, c = (//(1926o), h = l/(86 ), etc. This is the same as ignoring A S W in 
(|3~2^) and solving the resulting Hamilton- Jacobi equation as a local expansion. This local 
expansion is possible because the full solution to the Hamilton- Jacobi equation is just the 
Euclidean action evaluated on-shell and the classical theory is massive. Corrections to 
the coefficents can now be obtained iteratively by substituting the leading order solutions 
into A S W. 

For theories such as Yang-Mills which are classically massless the classical action does 
not have a local expansion, so that the semi-classical approach to solving ( |3"2"D is not 
viable. However the full theory is believed to be massive, |TT] , so that our expansions 



make sense, but we need a different approach to solving the resulting algebraic equations. 
We now propose a method in which the expansions in A are truncated at a certain order 
and the coefficents associated with the higher orders estimated using (]27f) . Consider, for 
example, the (^-independent contribution to A S W. Using 

V\ f A , 8M 166 2 A 2 2b nV ^X n 



3 

7T 



i/cT^/^io^T^ 1 ^ 1 ^ 1 ^' 2 <35) 



7 



(In the derivation of this formula we dropped the terms due to the delta-functions in (20 ) 
which only contribute a term proportional to 1/ y/s to this expression. Since this has no 
effect on the terms of higher order in s we ignore it here too). If we expand (y/sE — iq)~ l 
for large s in powers of l/y/s then for the high order terms the integral over the spectrum 
is suppressed by a large power of 1/E so that the dominant contribution will come from 
the lowest values of E. These correspond to the lightest particle one-particle states for 
which 

dp 



d^EmwmE)? = -fp 2 \(o r \m\E}\ 2 , P = vw^. (36) 

For large n 

J -J^P 2 f( E ) = J dpp 2 e^ n+1 ^ p2+m ^f(E) « J dpp 2 e-^ p2 / m2 m- 2 ^f(m) (37) 
This, together with Stirling's formula gives for large n 

^(-Mf^rV ^))- < 3s > 

n 6 \nm z J \ \nJ J 

with £ independent of A and n and proportional to m\(0 r \tp(0)\E = m)\ 2 . This shows that 
at high orders fl35l) is an alternating series. The error in truncating it at order \ N+1 / 2 
is proportional to (eA/ Nm 2 ) N+3 ' 2 / N 3 . If we set A = Nfi 2 /e then this error goes to zero 
with increasing N provided the ^-independent mass-scale /i satisfies (i < m, and at the 
same time A increases as we increase N. Similar, but more complicated calculations show 
that the coeflicents of ip 2 , p> 4 , p>' 2 behave like 




so that again if we truncate these series at order \ N+1 / 2 the error will go to zero as iV is 
increased. Furthermore, we obtain asymptotic values of the coefficents for large n, e.g. 

/ 1 \ n b ( 1 \ n c / 1 \ n 

hn ~ 2 -372' C ™ ~ 2 dn ~ e « ~ 2 dnV2 ( 40 ) 

with constant b, c, d. This suggests the following scheme for solving our local expansion 
of the Schrodinger equation. First truncate the coefficents of each linearly independent 
local functional of cp at, say, order \ N+1 / 2 taking A = Nfx 2 /e. This gives an infinite set 
of equations each involving a finite number of unknowns, but with an infinite number of 
unknowns overall. These can be solved with judicious use of ( j39D to estimate the discarded 
coefficents taking care to solve those equations for which the truncation error is smallest. 
In this approach there is no expansion parameter, but the approximation consists of 
working with large, but finite, values of N and A rather than the infinite values that occur 
in (|32"D. For the purpose of illustration only we take the unrealistically small value N = 0. 
Equating to zero the coefficents of 1, p 2 , p' 2 , p 4 , p> 2 p>' 2 p>" 2 keeping only the terms up to 
order vA gives six equations in the unknowns £ , bo, bi, b 2 , cq, ci, /o, fx- Our asymptotic 
formulae for the coefficents associated with higher orders in A are only valid for large n 
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but using them here for small n (after reinstating the factors of (n+l)/n that for large n 
were replaced by 1) we obtain c\ = —3y/2co/m 2 , fi = — 15f /m 2 ,b2 = —2y/2bi/(3y/3m 2 ). 
This gives us six equations in the six unknowns £, bo, b±, Co, fo, m. Thus we obtain 

C0 = ^(^-i|), ^ = ±( l + ^-M) (41) 



im 



where 60 satisfies 




g b V\-7riM 2 -^) lb 2 (lQ-17V2) + 2m 2 + y^-) =0 (42) 



Note that when the interaction is switched off, i.e. when g = 0, we obtain the free solution 
from M 2 — 4&o/ \pn = 0, as we should. Clearly it would be necessary to resort to numerical 
techniques to make progress with these algebraic equations. 

We can extend this approach to the calculation of the particle spectrum. Take the 
one-particle state corresponding to the lightest particle at rest to have the form of a pre- 
factor P mutliplied by ^ . For slowly varying fields this pre-factor will be the integral of 
a local function 

p = j dx (y + aiv? 3 + a 2 ip 5 + .. + hipip' 2 + 6 2 yV 4 + .. + Civ?V 2 + ..) (43) 
The Schrodinger equation equation for this state reduces to an equation linear in P 

h L ds s ~ v * (5< A -» P + I dx ^j^) + mP ) w = (44 > 

Again this may be reduced to a finite number of equations for a finite number of un- 
known coefncents by truncation and using formulae for the asymptotic behaviour of those 
coefficents derived from the study of A S P for large s. 

To summarise: the logarithm of the vacuum functional reduces to a local functional 
for fields that vary slowly on the scale of the lightest mass in the theory, but this does not 
satisfy the obvious form of the Schrodinger equation because expanding in slowly varying 
fields does not commute with removing a short distance cut-off. Instead we constructed 
the appropriate form of the equation for the case of ip 4 theory by re-summing the cut-off 
dependence using the analyticity properties of the Laplacian applied to the vacuum. This 
Schrodinger equation reduces the eigenvalue problem for the Hamiltonian to an infinite set 
of algebraic equations. We indicated how these may be solved in the usual semi-classical 
approximation and suggested a new approximation scheme that would also be viable for 
theories that are classically massless but acquire mass quantum mechanically such as the 
O(N) sigma model and Yang-Mills theories. 

Finally, I would like to acknowledge useful conversations with Hugh Osborn concerning 
the Schrodinger representation in field theory, and a grant from the Nuffield Foundation. 
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